Section: New Results
Towards the development of FFT-accelerated flexible fitting methods
Participants : Alexandre Hoffmann, Valerie Perrier, Sergei Grudinin.
We studied a set of new methods for non-rigid molecular fitting. The problem can be formulated as follows : Let and be two molecular structures (e.g. proteins). We are given , the electron density of and , the average positions of the atoms of . Assuming we can generate an artificial electron density from , our problem is to find a transformation of the atoms that minimizes the distance between and .
In image processing, this problem is usually solved using the optimal transport theory, but this method assumes that both densities have the same norm, which is not necessarily the case for the fitting problem. To solve this problem, one instead starts by splitting into a rigid transformation (which is a combination of translation and rotation) and a flexible transformation . Two classes of methods have been developed to find :
-
the first one uses optimization techniques such as gradient descent,
-
the second one uses the Fast Fourier Transform (FFT) to compute the Cross Correlation Function (CCF) of and .
We have developed several algorithms based on the FFT to find and we have developed two algorithms for flexible molecular fitting that are based on convex and non-convex optimization and the trust region methods. Our tests demonstrate that while one method gives good results for small deformations, the other gives good results for bigger deformations.
We have been also improving the current NMA method (which is essentially a model reduction technique), that is used in other tools such as the flexible fitting to small angle scattering profiles. Finally, we started the development of a method for a harder fitting/docking problem in which only electron density would be known. The basic idea would be to find the -diffeomorphism that minimizes the distance between and .
We developed several stand-alone C++ libraries to solve some of our problems including: